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ABSTRACT 

Using the Hubble Space Telescope, we have obtained the first color-magnitude diagram (CMD) to reach the 
main-sequence turnoff of the Galactic globular cluster NGC 6388. From a comparison between the cluster 
CMD and 47 Tucanae's, we find that the bulk of the stars in these two clusters have nearly the same age and 
chemical composition. On the other hand, our results indicate that the blue horizontal branch and RR Lyrae 
components in NGC 6388 are intrinsically overluminous, which must be due to one or more, still undetermined, 
non-canonical second parameter(s) affecting a relatively minor fraction of the stars in NGC 6388. 
Subject headings: stars: horizontal-branch — stars: variables: other — globular clusters: individual (47 Tu- 
canae, NGC 6388, NGC 6441) — globular clusters: general 



1. EVITRODUCTION 

NGC 6388 and NGC 6441 are two of the most intrigu- 
ing Galactic globular clusters (GC's). The integrated-light 
study by Rich, Minniti, & Liebert ( 1993) revealed a strong 
far-UV flux for these metal-rich ([Fe/H] ~ -0.60 and -0.53, 
respectively;iHarris 1996) bulge GC's. The far-UV flux in re- 
sol vedGC^sJ^ominatedb^JioUiorizontal branch (HB) stars 
(e.g., lDorman."b'Connell. & Rood 119951) . especially when 
rare UV-bright post-asymptotic giant branch stars are not 
present. Accordingly, the most likely explanation for the far- 
UV flux in NGC 6388 and NGC 6441 was immediately rec- 
ognized to be hot HB stars. However, while the UV- upturn 
phenomenon in elliptical galaxies is often attributed to blue 
HB stars (e.g., Dorman et al. 1995), no resolved metal-rich 
GC had been known with a blue HB morphology. The ten- 
dency for metal-rich GC's to have red HB's while metal-poor 
GC's have predominantly blue HB's reflects the "first param- 
eter" of HB morphology. As a consequence, the presence of 
blue HB stars in NGC 6388 and NGC 6441 would represent 
an example of the so-called "second-parameter (2"''P) phe- 
nomenon." 
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The presence of blue HB stars extending almost as faint in 
V as the turnoff (T O) point in both NGC 6388 and NGC 6441 
was confirmed bv iRic h et al.' fl997), who presented Hubble 
Space Telesc ope (HST ) photometry for both these GC's from 
the survey bv lPiotto'et al.i(i2002 ). A remarkable feature of the 
published diagrams is the presence of a strongly sloped HB at 
colors where other GC's ha ve a much m ore nearly "horizon- 
tal" HB. As emphasized bv lSweigart & Catelan C1998. here- 
after SC98), such a sloped HB cannot be simply the result 
of an older age or of enhanced mass loss along the red giant 
branch (RGB): while these are able to move a star horizontally 
along the HB, neither is able to increase the luminosity of a 
blue HB star compared to the red HB or RR Lyrae stars. Like- 
wise, while strong differential reddening might explain the 
sloping HB of a red HB cluster, it obviously cannot produce 
RR Lyrae and blue HB stars. SC98 conclude therefore that 
non-canonical 2"'^P candidates must be at play in NGC 6388 
and NGC 6441. 

However, th ese conclusions were challenged by 

iRaimondo et al.l i2002l hereafter R02), who computed 
models with non-standard values of the chemical abundance 
and mixing length parameter. Some of their models did reveal 
sloped HB's, but only as a consequence of an anomalously 
faint red HB (in V), together with a blue HB having a V-band 
luminosity consistent with the canonical models (see, e.g., 
their Fig. 2). 

Additional insights are provided by stellar variabil- 
ity and spectroscopic studies. Silbermann et al. ( 1994|, 
Lavdenetal. (1999), Pritzl et al. (2000, 2001, 2002, 2003), 
and .Corwin et al.. C2006.) have shown that the RR Lyrae vari- 
able stars in these GC's, which occupy the normally "hori- 
zontal" part of the HB, have much longer periods than field 
RR Lyrae of similar metallicity, thus strongly suggesting that 
they are intrinsically more lu minou s (SC9 8). Moreover, the- 
oretical calculations by Pritzl et al.l (!2002l) have shown that, 
contrary to the suggestions by Ree et al. (2002), the RR Lyrae 
components in both clusters cannot be explained in terms of 
evolution away from a position on the blue HB — and neither 
can the sloping nature of the HB be reproduced in this way. 
On the other hand, the first spectroscopic measurements of the 
gravities of blue HB stars in both N GC 6388 and NGC 6441 
(iMoehler, Sweigart, & Catelan 1999) revealed surface gravi- 
ties that are higher than predicted by even the canonical mod- 
els, thus arguing against an anomalously bright blue HB -i- RR 
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Fig. 1 .— The first deep CMD of NGC 6388 clearly revealing the TO point 
of the cluster {red X symbols) is overplotted on the 47 Tuc CMD (plus signs) 
in the V, B—I plane. To produce this plot, the 47 Tuc CMD was shifted by 
+3.2 in V and by +0.85 in B—I, so as to ahgn their red HB components. 



Lyrae component in these clusters. However, a recent re- 
assessment of the spectroscopic gravities of blue HB stars in 
NGC 6388 by Moehler & Sweigart (2006) indicates that the 
actual gravities should, in fact, be lower than the canonical 
values. It appears that the 1999 values must have been in error 
by a substantial amount, probably due to unresolved blends 
in the crowded inner regions of these massive {My ^ -9.5; 
IHarris 1996) GC's. 

In an effort to shed light on this puzzling situation, we have 
made use of the data obtained for NGC 6388 in the course 
of our snapshot HST program to study stellar variability in 
the cluster, and also of archival data, to produce its deepest- 
ever color-magnitude diagram (CMD). In §2 we describe this 
dataset and reduction procedures. In §3 we compare our CMD 
with 47 Tucanae's (NGC 104). We close in §4 by discussing 
the implications of our results for our understanding of the 
origin of the peculiar HB morphology of NGC 6388. 

2. OBSERVATIONAL DATA AND REDUCTION PROCEDURES 

The NGC 6388 data used in this paper were obtained un- 
der program SNAP-9821 (PI B. J. Pritzl), which used 
the Wide-Field Channel of the Advanced Camera for Surveys 
(ACS) to obtain six B (F435W), V (F555W), / (F814W) ex- 
posure triptychs on separate dates ranging from 2003 October 
to 2004 June. In addition, we have retrieved data from the 
HST Archives, as obtained under GO-9835 (PI G. Drukier). 
These consist of 12 V and 17 / exposures obtained with the 
High-Resolution Channel of ACS on 2003 October 30. 

The results employed here for 47 Tuc were obtained from 
the ground-based imagery used by Stetson to define his sec- 
ondary photometric standards. They consist of (125, 136, 94) 
images in (B, V, /) from 12 distinct observing runs; these pho- 
tometric indices should be on the system of Landolt ( 1992) 
to well under 0.01 mag. Stetson's ground-based data also in- 
clude (84, 137, 99) images in (B, V, I) for NGC 6388. Al- 
though these data are not included in the plots in this paper. 
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Fig. 2. — Similar to FigureQ but using the reddening-independent quantity 
Q instead of V, enforcing a match of the TO colors of the two clusters, and 
zooming in around the TO/SGB level. 



they were used to establish a network of photometric stan- 
dards in the cluster field which could be used to establish ac- 
curate photometric zero-points for the ACS images. Color 
transformations for the ACS data were based on these lo- 
cal standards and similar ACS-ground-based comparisons for 
47 Tuc, NGC 2419, NGC 6341 (M92), and NGC 6752. 

The data were reduced in the standard manner, using the 
DAOPHOT-ALLFRAME software packages, following com- 
monly understood reduction procedures (e.g.. Stetson 198"^ 
1990, 1994). A complete description of our dataset, reduction 
and calibration procedures will be described in a future paper 
(Stetson et al. 2006, in preparation). 

Before closing, we note that the ACS filter bandpasses are 
not id entical to the standard B, V, I bandpasses of LandoQ 
( 1992) — but then again, even among the various ground- 
based observing runs where the clusters were observed, the 
filters and detectors are not identical: bandpass mismatch is 
an unavoidable complication in filter photometry when one 
does not have a privately owned photometer. Even with a pri- 
vate photometer, bandpasses can drift with atmospheric con- 
ditions, or as filters, detectors, and mirror coatings age. We 
have removed the effects of mismatch to second order by in- 
cluding linear and quadratic color terms, both in the transfor- 
mation equations that relate our ground-based observations to 
Landolt's photometric system and in the equations that relate 
the ACS magnitudes to our ground-based system. Extensive 
experience suggests that the residual effects of bandpass mis- 
match in these filters reach an irreducible minimum scatter of 
0.01 - 0.02 mag on a star-by-star basis, and considerably less 
than this on an ensemble average of many stars. 

3. THE FIRST DEEP CMD OF NGC 6388 

Our deep V, B-I CMD of NGC 6388 is shown in Figure[T] 
(red X symbols), overplotted on the 47 Tuc CMD (plus signs). 
To obtain this plot, we have applied shifts of -1-3.2 in V and 
-1-0.85 in B-I to the 47 Tuc data, in order to register its red HB 
to NGC 6388's. The magnitude and color of the TO in each 
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cluster were determined by an iterative robust numerical fit of 
a parabola to the data in a magnitude range ±0.5 mag of the 
TO; we obtained (V, B-/)to = (20.93, 2.04) for NGC 6388, 
and (17.70, 1.27) for 47 Tuc. Note that we obtained clean, 
sharp, well-populated sequences in the CMD by plotting only 
the stars with the best photometry. Our selection criteria wiU 
be described in Stetson et al. (2006). 

3.1. The Main Sequence TO and Sub giant Branch (SGB) 

The TO points and SGB's of 47 Tuc and NGC 6388 coin- 
cide in brightness rather well once their red HB's are regis- 
tered. This suggests that any age difference between the two 
clusters is at the level of Af < 1 Gyr, assuming that their red 
H B's have the same absolute l uminosity. 

IStetson. VandenBere. & Bolte ( 1996) and 

[VandenBerg. Larson. & De Pr opris ( 1998J computed 
isochrones for different values of [Fe/H], [a/Fe], he- 
lium abundance Y, the mixing length parameter Omlt , and 
rotation rates uj. As is already well known, the difference in 
color between the base of the RGB and the main-sequence 
TO is a good relative age indicator, due to its small depen- 
dence on [Fe/H], [a/Fe], and uj. On the other hand, these 
calculations also show that this color difference, as well as 
the detailed shape of the SGB, do present some dependence 
on Y and omlt- 

To perform a more meaningful differential comparison that 
minimizes the effects of differential reddening in NGC 6388, 
we have replaced V with the reddening-independent quan- 
tity Q = V- 1.153 (B-/), and then registered the TO points 
of the two clusters in color Although numerical fitting of 
parabolas to 2 as a function of B-I at the TO's of the two 
GC's suggests that their bluest colors differ by 0.77 mag, we 
still performed a sanity check by directly overplotting the two 
CMDs with a range of horizontal shifts, and concluded that 
A(B-I) ^ 0.80 ±0.05 (1(7 error bar) appeared best to the 
eye (see Fig. |2j. We should note that the detection limit in 
these data is nea r Q « 19.5, and that past experience (e.g., 
IStetson & Harrisiil98 8) leads us to expect that the observed 
colors may be biased too blue by ~ 0. 1 mag as this limit is 
approached. Accordingly, we have been careful to restrict 
our comparisons to the immediate vicinity of the TO itself, 
Q ~ 18.2, where any systematic bias is expected to be much 
smaller. 

Figure |2] shows that, when the TO points of the clusters 
are registered, the difference in color A(Z?-/)!^§^ between the 
base of the RGB and the TO is very similar for both clusters. 
This quantity is only slightly larger for NGC 6388 than for 
47 Tuc, at a level ~ +0.010 ±0.035 mag (Ict error bar). If 
interpreted in terms of an age difference, this translates into 
Af ~ -0.5 ± 1.6 Gyr (NGC 6388 being younger). Assum- 
ing that the clusters have similar ages, [Fe/H], and [a/Fe], 
which appears most consistent with the TO/red HB ("vertical 
method") and the RGB data (see also the next subsection), one 
also finds AY ~ -0.014 ± 0.024 (in the sense that NGC 6388 
should have a lower Y), or AaMLT — -0.07 ± 0.12 (again in 
the sense that NGC 6388 should have a lower aMLx)- Such 
a difference in Y between the two GC's would be the oppo- 
site of that needed to explain the NGC 6388 HB morphology. 
Stronger constraints on Y and aMLX variations between the 
two clusters are derived from their RGB properties (see §3.2). 

3.2. The Red Giant Branch 

The RGB's of 47 Tuc and NGC 6388 have very similar 
morphologies, according to Figure ^ The larger scatter in 



the NGC 6388 CMD than in 47 Tuc's might be due in part to 
photometric errors and in part to differential reddening, thus 
making it difficult to determine if the stars scattered toward 
the red of the main RGB in NGC 6388 represent a minor 
metal-rich component. A sizeable metal-poor component is 
clearly not present (see also R02). A small deviation of the 
bulk of the brighter NGC 6388 RGB stars towards redder col- 
ors compared with the 47 Tuc CMD, if real, would suggest 
that NGC 6388 might be slightly more metal -rich (i.e., by 
< 0.25 dex) than 47 Tuc ([Fe/H] = -0.76 dex; 'Harris 1996^ 
When their red HB's are registered, one also finds that their 
TO colors become somewhat offset (with NGC 6388's be- 
ing bluer) and that their RGB's actually cross just above the 
HB level (Fig.Q. Theoretical isochrones show that these ef- 
fects are to be expected if NGC 6388 is more metal-rich than 
47 Tuc at the level indicated above. 

The "bump" in the RGB luminosity function lies at a V w 
17.70 in NGC 6388, and 14.50 in 47 Tuc. Thus, the TO- 
to-bump V-magnitude differences are 3.26 and 3.20, respec- 
tively. If real, this difference would suggest that NGC 6388 
is more metal-poor (by < 0.1 dex) than 47 Tuc. On the other 
hand, the values of AV^^™ measured by Zoccali et al. ( 199^ 
differ by only 0.03 mag between 47 Tuc and NGC 6388; this 
could be produced by a very small metallicity difference, by 
a Ay ~ 0.007 (NGC 6388 being more helium-rich; see Fig. 5 
in R02), and/or by a small difference in age of Af ~ 1 Gyr 
(NGC 6388 being younger). 

Note, finally, that near-infrared photometry (e.g., 
Froseletal. 2001; Valenti, Ferraro, & Orialia 2004) shows 
that NGC 6388 and NGC 6441 have, if anything, slightly 
bluer RGB's than 47 Tuc, contrary to what would be expected 
in the R02 scenario. Indeed, Ferraro et al. find a normal value 
of aMLT « 2.0 for both NGC 6441 and 47 Tuc (NGC 6388 
is not in their studied sample). The lack of a large differ- 
ence in aMLT and [Fe/H] between the two GC's, while not 
unexpected, already rules out the R02 scenario, in which a 
combination of high metallicity and low aMLT is invoked, 
leading to an underluminous red HB as mentioned in § 1 . 

3.3. The Red Horizontal Branch 

An overluminous red HB, as would be implied by a 
high primordial Y (SC98), is ruled out if, as appears 
likely, the two clusters have similar ages and metallici- 
ties. The HB luminosity function is also known to be 
affected by the stars' evolutionary parameters; in par- 
ticular, a higher Y in NGC 6388 should produce more 
luminosity evolution away fr om the zero-age HB (e.g., 
iDorman. VandenBerg. & Laskarides ,1989.) . which is clearly 
not present in the observed CMD. Likewise, the fact that 
there is no significant component fainter than the bulk of the 
NGC 6388 red HB stars suggests that any metal-rich compo- 
nent in this cluster (i.e., with [Fe/H] > -0.5 dex) should be 
minor. 

Figure |3l shows a more detailed comparison of the CMDs 
of NGC 6388 and 47 Tuc around the red HB region, again 
using the reddening-independent quantity Q. The 47 Tuc Q 
values were shifted by AQ = +2.22 mag. Adding random 
scatter at a level of cr(B-I) ~ 0.03 mag to the 47 Tuc CMD 
allows one to reproduce the overall appearance of the RGB 
and the red HB quite well, thus implying a AF < 0.03 (in 
the sense of NGC 6388 having a higher Y) based on the lu- 
minosity width of the HB. The main difference between the 
two HB's is the tendency for the bluer of the NGC 6388 red 
HB stars to be slightly brighter than the average line that de- 
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Fig. 3. — Same as in Figure^ but using Q instead of V and focusing around 
the red HB region. 



fines the 47 Tuc red HB. This, along with the fact that the TO 
points and SGB's of the two GC's match well in brightness 
when the two red HB's are registered, strongly suggests that 
both the RR Lyrae and blue HB components (and, to a lesser 
extent, also the bluer red HB stars) of NGC 6388 are indeed 
intrinsically overluminous with respect to the 47 Tuc red HB. 

3.4. The Asymptotic Giant Branch "Clump" 

Well-populated AGB "clumps" are seen i n both the 47 Tuc 
and NGC 6388 CMD's. 'Ferrar o et alJ ( fT999l) have shown that, 
the bluer the HB morphology of a GC, the less pronounced 
the resulting AGB clump. Therefore, the majority of the stars 
in this phase originate from red HB stars. It is interesting 
to note that the difference in V magnitude between the AGB 
clump and the red HB is basically indistinguishable between 
NGC 6388 and 47 Tuc. Unfortunately, this does not provide 



us with strong constraints on the difference in Y or metallicity 
between the two GC's, since the difference in magnitude be- 
tween the AGB clump and the red HB is not very sensitive to 
these evolutionary parameters (e.g.. Bono et al. 1995) . 

4. DISCUSSION 

In the present Letter, we have shown that, apart from the 
blue HB and RR Lyrae components, the CMD's of NGC 6388 
and 47 Tuc are very similar, thus strongly suggesting that the 
bulk of the stars in these two clusters are very similar Dif- 
ferences in age, metallicity, [a/Fe], Y, and Q!mlt between the 
two GC's, if present, should be small. Our results suggest 
that the red HB component of NGC 6388 is neither underlu- 
minous (as in the "canonical tilt" scenario of R02) nor over- 
luminous (as in the high-F scenario of SC98), except for its 
bluest stars. On the other hand, both the RR Lyrae and blue 
HB stars in NGC 6388 are significantly overluminous com- 
pared to field RR Lyrae stars of similar [Fe/H], thus explain- 
ing NGC 6388's anomalously long RR Lyrae periods. The 
lack of a sizeable luminosity difference between the red HB's 
of NGC 6388 and 47 Tuc indicates that the 2"''P that leads to 
the production of an overluminous RR Lyrae + blue HB com- 
ponent in NGC 6388 must be non-canonical in nature — that 
is, it must be neither age nor RGB mass loss, or else a sloped 
HB would not result (SC98). However, since only a relatively 
small fraction (~ 17%; Catelan 20 06) of the HB stars are in 
the RR Lyrae strip or on the blue HB, only a minor fraction 
of the cluster stars should be affected by this non-canonical 
2""*?, which could be either Y or/and the helium-core mass at 
the He-flash. (A fraction of the red HB stars, especially the 
bluer 25% or so, could also be affected, though to a much 
smaller extent.) Additional, detailed studies of the blue HB, 
RR Lyrae, and main sequence components will be required 
before we are in a positio n to conclusiv ely decide what 2"''P(s) 
is (are) responsible (see lCatelanl2006l for a recent review). 
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